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(Dated: February 15, 2008) 

Abstract 

We exploit the quantum coherence between pair-produced D° and D° in -0(3770) decays to study 
charm mixing, which is characterized by the parameters x and y, and to make a first determination 
of the relative strong phase 5 between doubly Cabibbo-suppressed D° — ► K + ir~ and Cabibbo- 
favored D° — > i"C + 7r~. We analyze a sample of 1.0 million D°D° pairs from 281 pb _1 of e + e~ 
collision data collected with the CLEO-c detector at E cm = 3.77 GeV. By combining CLEO-c 
measurements with branching fraction input and time- integrated measurements of Rm = {x 2 +y 2 )/2 
and i?ws = r(D° — > i^ + 7r _ )/r(l) — > K + 7r~) from other experiments, we find cos (5 = 1.03^gi7 ^ 
0.06, where the uncertainties are statistical and systematic, respectively. In addition, by further 
including external measurements of charm mixing parameters, we obtain an alternate measurement 
of cos 5 = 1.10 ± 0.35 ± 0.07, as well as x sin 5 = (4.4+^ ± 2 .9) x 10~ 3 and 5 = (22+^ ± 9 n )°. 

PACS numbers: 12.15.Ff,13.20.Fc,13.25.Ft,14.40.Lb 
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I. INTRODUCTION 



In the Standard Model, P°-P° mixing is suppressed both by the Glashow-Iliopoulos- 
Maiani mechanism [1] and by Cabibbo-Kobayashi-Maskawa matrix elements [2], although 
sizeable mixing could arise from new physics [3]. Charm mixing is conventionally described 
by two small dimensionless parameters: 

M 2 - Mi 

x = 2-== =-i 1 

r 2 + r\ 1 ; 

1/ = p — p , (2) 

1 2 + 1 1 

where M 12 and r 12 are the masses and widths, respectively, of the neutral P meson CP 
eigenstates, Pi (CP-odd) and P 2 (CP-even), which are defined as follows: 

\D°) + IP ) 
IP ) - |P°) 

W ^ (4) 

assuming CP conservation. The mixing probability is then denoted by R M = (a; 2 + 2/ 2 )/2, 
and the width of the P° and P° flavor eigenstates is T = (I\ + r 2 )/2. 

By focusing on P° decay times, recent experiments have made precise measurements 
of P°-P° mixing parameters [4-8] that highlight the need for information on the relative 
phase between the Cabibbo-favored decay P° — > K~n + and the doubly Cabibbo-suppressed 
(DCS) decay P° — > PJ~7r + . Direct measurements of y come from comparing decay times in 
P° — > K~7i + to those in P° transitions to the CP-even eigenstates K + K~ and 7r + 7r~. Time- 
dependent studies of the resonant substructure in P° — > PJ^7t + 7t~ provide x as well as y. 
In contrast, an indirect measure of y is provided by the "wrong-sign" process P° — > P + 7r~, 
where interference between the DCS amplitude and the mixing amplitude manifests itself in 
the P° decay time distributions. These analyses are sensitive to y' = y cos 5 — x sin S, where 
—5 is the relative phase between the DCS amplitude and the corresponding Cabibbo-favored 
P° — > P + 7r~ amplitude: (P + 7r~|P°)/(P' + 7r~|P ) = re~ lS . We adopt a convention in which 
5 corresponds to a strong phase, which vanishes in the SU(3) limit [9]. The magnitude r of 
the amplitude ratio is approximately 0.06. 

Measurements of y and y' have both attained a precision of 0(1O~ 3 ). However, because 5 
has not previously been measured, these separate determinations of y and y' have not been 
directly comparable. Thus, even a modest measurement of 5 can significantly improve the 
overall knowledge of charm mixing parameters. 

In this article, we present a first determination of 5 that takes advantage of the correlated 
production of P° and P° mesons in e + e~ collisions produced at the Cornell Electron Storage 
Ring and collected with the CLEO-c detector. If a collision produces no accompanying 
particles, the P°P° pair is in a quantum-coherent C = — 1 state. Because the initial state 
(the virtual photon) has J PC = 1 , there follows a set of selection rules for the decays of 
the P° and P° [3, 9-18]. For example, both P° and P° cannot decay to CP eigenstates with 
the same eigenvalue. On the other hand, decays to CP eigenstates of opposite eigenvalue are 
enhanced by a factor of two. More generally, final states that can be reached by both P° and 
P° (such as K~tc + ) are subject to similar interference effects. As a result, the effective P° 
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branching fractions in this D°D° system differ from those measured in isolated D° mesons. 
Moreover, using time-independent rate measurements, it becomes possible to probe D°-D° 
mixing as well as the relative strong phases between D° and D° decay amplitudes to any 
given final state. 

We implement the method described in Ref. [19] for measuring y and 5 using quantum 
correlations at the ?/>(3770) resonance. Our experimental technique is an extension of the 
double tagging method previously used to determine absolute hadronic P-meson branching 
fractions at CLEO-c [20]. This method combines yields of fully reconstructed single tags 
(ST), which are individually reconstructed D° or D° candidates, with yields of double tags 
(DT), which are events where both D° and D° are reconstructed, to give absolute branching 
fractions without needing to know the luminosity or D°D° production cross section. Given 
a set of measured yields, efficiencies, and background estimates, a least-squares fitter [21] 
extracts the number of D°D° pairs produced (A/") and the branching fractions (£>) of the 
reconstructed D° final states, while accounting for all statistical and systematic uncertainties 
and their correlations. We employ a modified version of this fitter that also determines y, 
x 2 , r 2 , r cos 5, and rx sin 8 using the following categories of reconstructed final states: K~n + 
(/), K + ix~ (/), CP-even (S+) and CP-odd (S-) eigenstates, and semileptonic decays (e 1 * 1 ). 
With CLEO-c measurements alone, r 2 is not determined with sufficient precision to extract 5. 
Therefore, we also incorporate measurements of branching fractions and mixing parameters 
from other CLEO-c analyses or from external sources. We neglect CP violation in D decays, 
which would entail a slight correction to the mixing signal. 

The paper is organized as follows. In Section II, we review the formalism of quantum- 
correlated D°D° decay. Section III describes the event selection criteria and D reconstruction 
procedures. The external measurements used in the fit are summarized in Section IV. 
Systematic uncertainties, which are also input to the fit, are discussed in Section V. Finally, 
we present and discuss our main fit results in Section VI. In the appendix, we provide 
information for use by other experiments. 

II. FORMALISM 

To first order in x and y, the C-odd rate T D o£)o(i, j) for D°D° decay to final state {i, j} 
follows from the anti-symmetric amplitude M.if 

T D05 o(i,j) oc M% = \A t A 3 - AiAj\ 2 

= mD^ijlD^-i^ijlD^ 2 , (5) 

where Ai = (i\D°) and A\ = (i\D°). These amplitudes are normalized such that B K - n + xs 
A 2 r _ ?r+ (1 +ry cos 5 + rx sin 5), B s± ~ A 2 s± (l^y), and B e ps A 2 e . The total rate, T^o^o, is the 
same as for uncorrelated decay, as are ST rates. However, unlike the case of uncorrelated 
P°P°, we can consider the C-odd D°D° system as a P1P2 pair. If only flavored final states 
are considered, as in Ref. [20], then the effects of quantum correlations are negligible. In 
this analysis, we also include CP eigenstates, which brings additional sensitivity to y and 
5, as demonstrated below. 

Quantum-correlated semileptonic rates probe y because the decay width does not depend 
on the CP eigenvalue of the parent D meson, as this weak decay is only sensitive to flavor 
content. However, the total width of the parent meson does depend on its CP eigenvalue: 

Ti = r(l =p y), so the semileptonic branching fraction for D x or D 2 is modified by 1 ± y. 

2 
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If we reconstruct a semileptonic decay in the same event as a D2 — > S+ decay, then the 
semileptonic D must be a Di. Therefore, the effective quantum-correlated D°D° branching 
fractions (jF cor ) for {S±,e} final states depend on y: 



-rcor ^ 



2B 5± B e (l±y), 



(6) 



where the factor of 2 arises from the sum of e + and e~ rates. Combined with estimates 
of B e and B s± from ST yields, external sources, and flavor-tagged semileptonic yields, this 
equation allows y to be determined. 

Similarly, if we reconstruct aD^ K~n + decay in the same event as a D2 — > S+, then 
the K~7i + was produced from a Di. The effective branching fraction for this DT process is 
therefore 



-77Cor 

J ~S+,K--k j < 



KS+lD^iK-n+lDj] 2 

Bs+B K - n +(l + y)(l - ry cos 5 - r:rsin5)|l + re~ i5 \ 2 
Bs + Bk-k+(1 + — rycosS — rxsin<5)(l + 2r cos5 + r 2 ) 
Bs + B K -n+ (1 + 2r cos 5 + R ws + y) , 



(7) 



where i?ws is the wrong-sign rate ratio, which depends on x and y because of the interference 
between DCS and mixing transitions: i?ws = r(D° — > fT"7r + )/r(D° — > K~tt + ) = r 2 + ry f + 
i?M- In an analogous fashion, we find ^ r 5° r ^-- 7r + ~ Bs_Bk--k+(^ + Rws — 2rcos5 — y). When 
combined with knowledge of Bs + , y, and r, the asymmetry between these two DT yields 
gives cos 5. In the absence of quantum correlations, the effective branching fractions above 
would be B s± B K - n +(l + Rws)- 

More concretely, we evaluate Eq. (5) with the above definitions of r and 8 to produce the 
expressions in Table I. In doing so, we use the fact that inclusive ST rates are given by the 
incoherent branching fractions since each event contains one D° and one D°. Comparison of 
jF cor with the uncorrected effective branching fractions, jF unc , also given in Table I, allows 
us to extract r 2 , rcos<5, y, x 2 , and rxsin5. Information on Bi is obtained from ST yields 
at the -0(3770) and from external measurements using incoherently produced D° mesons. 
These two estimates of Bi are averaged by the fitter to obtain jF unc . 

Using only ST and DT yields at the -0(3770), we can determine y and cos S from the 
following double ratios, obtained by manipulating the expressions in Table I: 



N 



fe 



AN 



N s . 
N s _ 



N. 



s+ 



2r cos S + y 



1 

4 
N 



f 



N s+e N s _ 
Ns. e N s+ 



N s+ 

Ns^eNs, 

' N s+e N s _ 



ff 



AN 1 



N, 



N 



s+ 



Ns.f 



'-N f N 



f^S+S- 



Ns +f 
1 



N f s-N s+ N fS+ N s _ 



(8) 
(9) 
(10) 

(11) 



where N denotes an efficiency-corrected background-subtracted ST or DT yield. Note that 
semileptonic yields are essential for separating y and cos 5. Including external branching 
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TABLE I: Correlated (C-odd) and uncorrelated effective D°D° branching fractions, J- COT and J 7 ™ 110 , 
to leading order in x, y and Rws, divided by B{ for ST modes i (first section) and B{Bj for DT 
modes {i,j} (second section). Charge conjugate modes are implied. 



Mode 


Correlated 


Uncorrelated 


K-1T+ 


l + Rws 


l + Rws 


S + 


2 


2 


S- 


2 


2 




i?M 




K-TT + ,K + 7T- 


(1 + Rws) 2 — 4r cos 5{r cos 5 + y) 


l + R-ws 


K-TT+,S + 


1 + i? ws + 2r cos 5 + y 


l + Rws 




1 + Rws — 2r cos 5 — y 


l + Rws 


K~ir + , e~ 


1 — ry cos 5 — rx sin 5 


l 


s+, s+ 





1 


S-,S- 





1 


S + ,S- 


4 


2 


S + ,e~ 


i + y 


1 




i-y 


1 



fractions provides additional ways to determine y and cos 5: 

2AfB s± B e 



y 



l - 



S±e 



2rcos5 + y = (1 + -Rws) 



N s+f /B s+ - N s _ f /Bc 
N s+f /B s+ + N s _ f /B< 



(12) 
(13) 



Although we neglect x 2 and y 2 terms in general, we report a result for x 2 as determined 
solely from Nff/Nfj. 



III. EVENT SELECTION AND RECONSTRUCTION 



Our current analysis uses 281 pb _1 of e + e~ — > -^(3770) data collected with the 
CLEO-c [22-26] detector. We also make use of a large Monte Carlo simulated sample 
of uncorrelated D°D° decays with an effective luminosity 40 times that of our data sam- 
ple, from which we estimate signal efficiencies, background contributions, and probabilities 
for misreconstructing a produced signal decay in a different signal mode (crossfeed). In 
these samples, we reconstruct the final states shown in Table II, with 7r° — > 77, r\ — > 77, 
Kg — > 7r + 7r~, and uo — > 7r + 7r~7r°. Because most K\ mesons and neutrinos are not detected, 
we do not reconstruct K^n and semileptonic ST modes; they are only included in DT 
modes, paired with a fully reconstructed D candidate. Below, we denote by S' + the subset 
of S + modes that are fully reconstructed: K + K~, 7r + 7r _ , and KgTT°TT°. For CP eigenstates, 
we choose modes with unambiguous CP content. In addition to two-body decays, we also 
include fT^7r 7r°, which is a pure CP-even eigenstate because the two identical 7r°'s must 
have even angular momentum in order to satisfy Bose symmetry. We neglect CP violation 
in K° decays. 
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TABLE II: D final states reconstructed in this analysis. 



Type 



Final States 



Flavored 




K-TT+, K+ir- 
K+K~, vr+vr-, K%ir°ir°, 

Inclusive Xe + v e , Xe~v e 



.0 



Standard CLEO-c selection critera for 7r ± , K ± , ir°, and K s candidates are described in 
Ref. [27]. In addition, for K s candidates, we impose \M(n + n^) — M K a \ < 7.5 MeV/c 2 , and 
we require the decay vertex to be separated from the interaction region with a significance 
greater than two standard deviations (<r). For the K s tt°tt° mode, we reject K s daughter 
tracks with ionization energy loss consistent with being kaons at the level of 2a in order to 
suppress misreconstructed D — > K~it + it Q decays, where the kaon is taken to be a charged 
pion, and a soft combinatoric n° candidate is incorporated into the D candidate. We accept 
lu candidates with |M(7r + 7r~7r°) — M w | < 20 MeV/c 2 . Reconstruction of rj — > 77 proceeds 
analogously to 7r° — > 77. In addition, we require |M(77) - M v \ < 42 MeV/c 2 . All invariant 
mass requirements correspond to approximately 3<r consistency with the nominal masses 
from Ref. [28]. 

1. Single Tags 

Reconstruction of all modes in this analysis, including DT modes, begins with fully 
reconstructed ST D candidates, which are identified using two kinematic variables that 
express momentum and energy conservation: the beam-constrained candidate mass M and 
the energy difference AE. These variables are defined to be 



where p^> and E D are the total momentum and energy of the D candidate, and E is the 
beam energy. Correctly reconstructed D candidates produce a peak in M at the D mass 
and in AE at zero. We determine ST yields by fitting the M distribution with the mode- 
dependent requirements on AE listed in Table III, which are applied to both ST and DT D 
candidates. The limits are set at approximately three standard deviations. Modes with 7T° 
or 77, which decay to two photons, have asymmetric limits to allow for partially contained 
showers in the electromagnetic calorimeter. 

For K ± 7i Zf , K + K~ , and 7r + 7r - ST modes, in events containing only two tracks, we suppress 
cosmic muons and Bhabhas scattering events by vetoing tracks that are identified as muons 
or electrons and by requiring at least one electromagnetic shower in the calorimeter above 50 
MeV not associated with the signal tracks. For K + K~ ST candidates, additional geometric 
requirements are needed to remove doubly radiative Bhabhas followed by pair conversion of 
a radiated photon. Also, we accept only one candidate per mode per event; when multiple 
candidates are present, we choose the one with smallest [Ai?]. 

The resultant M distributions, shown in Fig. 1, are fitted to a signal shape derived from 
simulated signal events and to a background ARGUS function [29]. The simulated signal 




(14) 
(15) 



AE = Ed — Eq, 
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TABLE III: Requirements on AE for D candidates. 



Mode 


Requirement (GeV) 


K-TT+ 


\AE\ < 0.0294 


K + TT 


\AE\ < 0.0294 


K+K- 


\AE\ < 0.0200 


7T + 7T~ 


\AE\ < 0.0300 




-0.0550 < AE < 0.0450 


K° S V 


-0.0710 < AE < 0.0450 


-0.0550 < AE < 0.0350 




|A£7| < 0.0250 



shape is both shifted and convoluted with a Gaussian smearing function to account for im- 
perfect modeling of the detector resolution and beam energy calibration. The width of the 
smearing function is allowed to float in the fit. The measured ST yields and efficiencies 
are given in Table IV. The yield uncertainties are statistical and uncorrected systematic, 
respectively. The latter arise from modeling of multiple candidates in simulation and vari- 
ations in the signal lineshape. Correlated systematic uncertainties are discussed separately 
in Section V. 

TABLE IV: ST yields and efficiencies input to the data fit. Yield uncertainties are statistical and 
uncorrelated systematic, respectively, and efficiency uncertainties are statistical only. 



Mode 


Yield 


Efficiency (%) 


K-TT+ 


25374 ± 166 ± 26 


64.70 ± 0.04 


K+ir- 


25842 ± 167 ± 26 


65.62 ±0.04 


K+K- 


4740 ± 71 ± 5 


57.25 ± 0.09 




2098 ± 59 ± 9 


72.92 ±0.13 


K° S TT°TT 


2435 ± 72 ± 16 


12.50 ±0.06 


K° s n° 


7523 ± 91 ± 17 


29.73 ± 0.05 


K° s r] 


1051 ±39 ±17 


10.34 ±0.06 




3239 ± 63 ± 7 


12.48 ± 0.04 



2. Fully Reconstructed Hadronic Double Tags 

Except for modes with K^ir , we form hadronic DTs by combining two ST candidates 
passing the above selection criteria. Multiple candidates are resolved after forming the DT 
candidates, not at the ST level. We choose one candidate per mode per event with M 
closest to the measured D° mass, where M = [M(D°) + M(D°)}/2. We extract signal 
yields by counting events in the two-dimensional M(D°) vs. M(D°) plane, as illustrated in 
Fig. 2. The signal region S is defined to be approximately three standard deviations in each 
dimension: 1.86 < M(D°) < 1.87 GeV/c 2 and 1.86 < M(D°) < 1.87 GeV/c 2 . Sidebands A 
and B contain candidates where either the D° or the D° is misreconstructed. Sidebands C 
and D contain candidates where both D° and D° are misreconstructed, either in a correlated 



8 



3950208-003 




1.83 1.84 1.85 1.86 1.87 1.88 1.89 1.83 1.84 1.85 1.86 1.87 1.88 1.89 



M (GeV/c 2 ) 

FIG. 1: ST M distributions and fits. Data are shown as points with error bars. The solid lines 
show the total fits, and the dashed lines show the background shapes. 

way (C) by assigning daughter particles to the wrong parent or in an uncorrelated way (D). 
Event counts in sidebands A, B, and C are projected into the signal region S using scale 
factors determined from integrating the background shape in the ST M fits. Contributions to 
sideband D are assumed to be uniformly distributed across the other regions. To account for 
systematic effects in the sideband definitions and in the extrapolation to the signal regions, 
we assign a 100% systematic uncertainty on the size of the sideband subtractions, which is 
much smaller than the statistical uncertainties in all cases. 

Table V gives the fully reconstructed DT yields and efficiencies input to the fit, and Fig. 3 
shows the corresponding M projections. Same-CP {S±, S±} modes are not included in the 
standard fit. 
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M (GeV/c 2 ) 

FIG. 2: Two-dimensional M distributions with signal (S) and sideband (A, B, C, D) regions 
depicted, for the sum of all {K^ir^, S' + } and {K^-k^, S-} modes listed in Table V. 



3. Semileptonic Double Tags 

Semileptonic DTs are partially reconstructed by combining a fully reconstructed hadronic 
ST with an electron candidate from the remainder of the event. The hadronic tags are 
selected using the same criteria as for STs, described in in Section III 1, with an additional 
mode-dependent M requirement listed in Table VI. The limits are set at approximately 
three standard deviations. Multiple tag candidates are resolved with minimal |A£?|, as for 
STs. All electron candidates in a given event are accepted. 

Electron candidate tracks are selected using the criteria described in Ref . [30] . Electrons 
are distinguished from hadrons via a multivariate discriminant that combines information 
from the ratio of the energy deposited in the calorimeter to the measured track momen- 
tum (E/p), ionization energy loss in the tracking chambers (dE/dx), and the ring-imaging 
Cerenkov counter (RICH). When paired with a f^ =F 7r =t tag, the kaon and electron charges 
must be the same. This requirement removes charge-symmetric backgrounds and cannot be 
used with unflavored S± tag modes. 

As illustrated in Fig. 4, the resulting electron momentum spectrum for each tag mode 
is fit to signal and background shapes fixed by the simulation; only the normalizations are 
allowed to float. We correct the simulated signal spectrum for relative bin-by-bin efficiency 
differences between data and simulation, which are measured with a high-statistics radiative 
Bhabha sample. The background components in the spectrum fit include misreconstructed 
tags, photon conversions and light hadron decays to electrons (mostly ir° Dalitz decays), weak 
decays in flight, and hadron misidentification. The latter component peaks at 700 MeV/c 
and mostly consists of K that escape the acceptance of the RICH; are indistinguishable 
from using dE/dx at this momentum. In the fit, the misreconstructed tag component is 
fixed from sidebands in M and AE. The other backgrounds contribute approximately 7% 
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FIG. 3: Sums of fully reconstructed DT M distributions, with charge conjugate modes implied. 
Data are shown as points with error bars. The shaded histograms show events outside the signal 
region. 



for K^7f ± tags and 20-30% for S± tags. 

Table VII gives the semileptonic DT yields and efficiencies input to the fit. The uncorre- 
lated systematic uncertainties are determined from yield excursions under variation of the 
signal and background shapes used in the spectrum fit. For the signal shape, we adjust the 
semileptonic form factor model used in the simulation. For the background shapes, we vary 
the allowed background composition by removing components with insignificant yields and 
also allowing for negative normalizations. 



4- Double Tags with K^ir 

DT modes with K^tt are reconstructed with the missing mass technique also used in 
Ref. [31]. A fully reconstructed ST candidate, selected as described in Section III 1, is 
combined with a tt° candidate, and we compute the squared recoil mass against the ST-7T 
system, M^ iss . Signal Kin decays produce a peak in M^ iss at M K o. Backgrounds from 
D — ► KgiT , 7r°7r°, and r/7r° are suppressed by vetoing events with additional unassigned 
charged particles, 77, or tt candidates. 

Figure 5 shows the resultant M^ iss in data. Signal yields are obtained from event counts 
in the signal region of 0.1 < M^ iss < 0.5 GeV 2 /c 4 . The contribution from combinatoric 
background is estimated from the 0.8 < M^ iss < 1.2 GeV 2 /c 4 sideband to be O(10 -3 ) and 
is subtracted. The residual background contributions from K^n and tjtt , which peak in 
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TABLE V: Fully reconstructed DT yields and efficiencies input to the data fit. Yield uncertainties 
are statistical and uncorrelated systematic (for sideband subtraction), respectively, and efficiency 
uncertainties are statistical only. 



Mode 




Yield 


Efficiency (%) 


K-tt h 


-,K-ir + 


2.0 ± 1.4 ±0 


36.1 ±3.4 


K-^ 


-,K + TT- 


600 ± 25 ± 5 


41.1 ±0.2 


K-*+ 


-,K+K~ 


71 ± 8 ± 1 


35.5 ±0.6 




", 7T + 7r~ 


24 ± 5 ± 1 


44.4 ± 1.1 




-,K° S TT°TT 


32 ± 6 ± 1 


8.0 ±0.3 


K-^ 




88 ± 9 ± 1 


18.4 ±0.3 


K-^ 


-,K° sV 


8.0 ±2.8 ±0.0 


6.0 ±0.3 


K-^ 




29 ± 5 ± 


8.7 ±0.2 


K + 7t- 


-,K+TT- 


2.0 ±1.4 ±0.0 


44.1 ± 3.3 


K + it- 


,K+K- 


54 ± 7 ± 


36.1 ±0.6 


K + it- 


", 7T + 7T~ 


25 ± 5 ± 1 


48.1 ± 1.1 


K + it- 


-,^7r°7r° 


33 ± 6 ± 


8.0 ±0.3 




-,K° S 7T° 

-,K° sV 


76 ± 9 ± 


18.6 ±0.3 


K+TT' 


9±3±0 


6.1 ±0.3 


K+TT- 




33 ± 6 ± 1 


8.0 ±0.2 


K+K 


- xg* 


39 ± 6 ± 1 


17.1 ± 0.6 


K+K 




7.0 ±2.7 ±0.0 


7.1 ±0.8 




20 ± 4 ± 


6.8 ±0.4 






13±4±0 


19.2 ± 1.2 


7T + 7r~ 


,K° s ri 


2.0 ±1.4 ±0.0 


8.1 ± 1.4 


7T + 7r~ 


,K° s uj 


7.0 ±2.7 ±0.0 


9.9 ±0.9 


K° s ir° 


ir°,K° s Tr 


14 ± 4 ± 


3.5 ±0.2 


K S*° 




4.0 ±2.0 ±0.0 


1.2 ±0.2 






4.0 ±2.0 ±0.0 


1.4 ±0.2 



the signal region, and from 7r°7r°, which peaks below the signal region, are estimated from 
simulated data and are subtracted by the fitter, as described in Section III 5. 

Table VIII lists the K^n DT yields and efficiencies input to the fit. There are no 
uncorrelated systematic uncertainties for these modes. 

5. Crossfeed and Peaking Backgrounds 

The yields in Tables IV, V, VII, and VIII include peaking backgrounds contributions and 
crossfeed among signal modes that are subtracted by the fitter. To do so, the fitter makes 
use of crossfeed probabilities and background efficiencies determined from simulated events, 
as well as branching fractions for peaking background processes [28] . These inputs are listed 
in Table IX. The correlated uncertainties among these contributions is also accounted for. 

In most cases, the peaking background processes produce the same final state particles 
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TABLE VI: Requirements on M for hadronic tags in semileptonic DT candidates 

Mode Requirement (GeV/c 2 ) 

K-7T+ 1.8585 < M < 1.8775 

K+TT- 1.8585 < M < 1.8775 

K + K~ 1.8585 < M < 1.8775 

7T+7T- 1.8585 < M < 1.8775 
K° s tt°tt 1.8530 < M < 1.8780 

K%ir° 1.8530 < M < 1.8780 

K^r] 1.8585 < M < 1.8775 

K^lo 1.8585 < M < 1.8775 



TABLE VII: Semileptonic DT yields and efficiencies input to the data fit. Yield uncertainties 
are statistical and uncorrelated systematic, respectively, and efficiency uncertainties are statistical 
only. 



Mode 






Yield 


Efficiency (%) 


X+e~ 




K-7T+ 


1128 ± 45 ± 13 


45.3 ±0.2 


X+e- 


Ve 


K-K+ 


128 ± 24 ± 5 


40.1 ±0.5 


X+e- 


Ve 


7T~7T + 


49 ± 12 ± 7 


50.3 ±0.9 


X+e- 


l>e 


^7r 7T 


37 ± 22 ± 12 


8.9 ±0.2 


X+e- 


l>e 


K° S TT° 


195 ± 24 ± 2 


21.4 ±0.3 


X + e- 


Ve 


K° S 7] 


28 ± 6 ± 4 


7.9 ±0.3 


X+e- 


Ve 


K°LJ 


50 ± 15 ± 6 


8.2 ±0.2 


X~e+ 


Ve 




1218 ± 47 ± 17 


45.9 ±0.2 


X-e+ 


Ve 


K-K+ 


102 ± 21 ± 10 


39.0 ±0.5 


X~e+ 


Ve 


TT~TT + 


40 ± 10 ± 4 


50.5 ±0.8 


X~e+ 


V e 




50 ± 15 ± 9 


9.6 ±0.2 


X~e+ 


Ve 


K° s n° 


189 ± 19 ± 7 


21.4 ±0.3 


X-e+ 


V e 


K° S V 


27±8±2 


7.7 ±0.3 


X~e+ 


V e 




49 ± 16 ± 4 


8.5 ±0.2 



as the signal modes to which they contribute. For these backgrounds, the DT contribution 
is assumed to occur at the same rate as for STs. Other backgrounds, indicated by asterisks 
in Table IX, arise from misreconstructed D decays. These backgrounds are present only 
in ST modes and do not contribute to DT modes because of the kinematic constraints of 
full event reconstruction. Background processes are identified using simulated D°D° events, 
where each D is uncorrelated with the other and decays generically. Background efficiencies 
are determined by observing the change in signal yield as each background contribution is 
removed from the simulated sample. 

We also adjust the peaking background estimates for quantum correlation effects. The 
background branching fractions and efficiencies mentioned above assume uncorrelated D 
decay, but the background estimate also depends on the type of decay of the other D. For 
instance, K^tc signals have crossfeed background from K^ir but only for Ktt and CP-even 
tags. In K®ir° candidates opposite CP-odd tags, K^ir cannot contribute because it is also 
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FIG. 4: Sums of electron momentum spectra for X e + u e and X + e v e . Data are shown as points 
with error bars, and the histograms show the signal and background components in the fits. 



CP-odd. For non-CP-eigenstate multi-body backgrounds, we assume equal CP-eigenstate 
and flavored content, with the CP-eigenstate content equally divided between CP-even and 
CP-odd. Variations in these assumptions give rise to systematic uncertainties, which are 
assessed in Section VI. In general, peaking backgrounds account for less than 1% of the 
measured yields, except in Kgcu modes (5-10%), K^n modes (1-2%), and ST Kgn n° 
(3%). 
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FIG. 5: Sums of M^ iss distributions for K^ir . The shaded histograms represent simulations of 
the peaking backgrounds D° — > 7r 7r°, K^tt , rjir , and K*°7r°. 



TABLE VIII: K^ir DT yields and efficiencies with statistical uncertainties input to the data fit. 

Mode Yield Efficiency (%) 

Klir°,K-Tr + 187 ± 14 34.0 ± 0.2 

K%ir°,K + ir- 179 ± 14 34.8 ± 0.2 

K°ir°, K°ir° 90 ± 10 15.5 ± 0.1 

K°Tr°,K°r] 8.6 ±3.1 5.4 ±0.1 
K<l7r ,K°uj 34 ±6 6.2 ±0.1 



IV. EXTERNAL MEASUREMENTS 

As discussed in Section II, external estimates of uncorrelated branching fractions can help 
determine y and 5. We include measurements from Refs. [28, 32] and from other CLEO-c 
analyses along with their full error matrix in the fit. For the CLEO-c analyses, we account 
for statistical correlations with the yields measured in this analysis. Tables X, XI, and XII 
shows the sources of these measurements, and there is no overlap among the measurements 
in these Tables. 



A. External Branching Ratio Measurements 

Because our precision on cos 5 is currently limited by our knowledge of the CP-eigenstate 
branching fractions, we include as many external measurements of these branching fractions 
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TABLE IX: Crossfeed probabilities and peaking background efficiencies and branching fractions 
input to the fit. Backgrounds marked by an asterisk (*) occur in only in STs, not DTs. 



Crossfeed 




Signal 




Jrrobabinty [/of 


K TV 




K TV 




U.Uyo ± U.UUo 


J\ TV 




K TV 




U.Uyz zb U.UUo 


TSO, , 
S 




KgTV TV 




n H/i c j_ n no i 
U.U4o zb U.Uol 






Kj^TV 




i no _i_ n no 
l.\3l ± U.Uo 


Background 


* 


signal 


"bkg {/o ) [28\ 


^mciency { /o) 


KgTV^TV 




tsO 

KgTV TV 


2.90 ± 0.19 


0.006 ± 0.012 


TV^TV 7T U 7T U 




KgTV TV 


u.9o ± u.uy 


0.029 ± 0.067 


/ - 4- /;k\ 
K TV^TV (*) 




tsO 

KgTV TV 


14.1 ± 0.5 


0.0012 ± 0.0024 


— > KgTV^TV (*) 




tsO 
KgTV TV 


7.22 ± 0.26 


0.048 ± 0.007 


+ — 


* 


T^O 

K s^ 


1.00 ± 0.0b 


0.071 ± 0.005 


„0„0 
p TV 




z^0„0 
KgTV u 


0.32 ± 0.04 


0.032 ± 0.033 


Generic JJ L> ( ) 


* 


KgTT 




0.00051 ± 0.00016 


Generic D^D {*) 




T^O 

K S^ 




0.00038 ± 0.00012 


K TV TV 






2.3 ± 0.5 


0.065 ± 0.002 


is*— + 
K TV^TV 




ISO 
KgLO 


1.00 ± 0.22 


0.135 ± 0.003 


i\ p 




S 




n A^^ -1- A AA1 


K*°p° 






1.50 ±0.33 


0.022 ± 0.004 


K^TV+TV'TV 




K° s uj 


1.1 ± 1.1 


0.405 ± 0.423 


7T + 7r~7r + 7T~7r 




K° s uj 


0.41 ±0.05 


0.006 ± 0.003 


K^TV + 




K° s uj 


1.12 ±0.31 


0.038 ± 0.005 


K?7T° 




K° s uj 


0.59 ±0.02 


0.030 ± 0.010 


rjTV 




K° l tv° 


0.056 ±0.014 


0.155 ±0.001 


TV°TV° 




K°TV° 


0.079 ± 0.008 


0.117 ±0.001 



as possible. Except for the inclusive B(KgTv°TV ) and B^K^tv ), all CP-eigenstates in this 
analysis have previous branching ratio measurements with respect to other reference modes. 

For KgTv , Kgi], and KgU, we make use of the global branching fraction fit performed in 
Ref. [32]. Because B(K~~tt + ) is a free parameter in this fit and in the current analysis, we also 
include Ref. [32] 's fit result for this mode in order to properly account for the correlations 
among these branching fractions. We do not use Ref. [28] because this compilation includes 
a CLEO-c measurement of £>(fT~7r + )[2u], which is correlated with the measurements in this 
analysis. 

Previous experiments have measured Rkk = B(K~~ K + )/B(K~tv + ) and R n7T = 
B{tv~tv + ) /B(K~tv + ) simultaneously [28], so these two quantities are correlated both sta- 
tistically [via the common denominator, B(K~tv + )] and systematically. It is common for 
these experiments to also report a value of Rkk/R^ m addition to Rkk and R nn sepa- 
rately. Using the dominant measurements of all three quantities [40, 48-50], we compute a 
weighted average correlation coefficient of 0.30 between the Rkk and R nn values quoted in 
Ref [28]. For R^, we remove the CLEO-c measurement [51] from the average because it 
was based on the same dataset as this analysis. 

The CLEO-c D — > K% l tv° analysis [31], provides additional information on B(K^tv°) 
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from the K~7i + 7i° and K~7i + 7i~7i + tag modes, which are not used in this analysis. Sys- 
tematic correlations between the K^tt yields in this analysis and these additional i3(i^°7r°) 
measurements is taken into account. We also include the statistical correlation between 
i3(-K"°7r°) and B(Kgir°); knowledge of the latter is required to correct for the effect of quan- 
tum correlations in the former. 

Table X summarizes the external branching fraction inputs to the fit. 

TABLE X: External branching ratio measurements and their correlations used in the data fit. An 
asterisk (*) indicates removal of overlapping CLEO-c measurements. 



Parameter Value (%) Correlation Coefficients 

B{K-ir+) 3.81 ± 0.09 [32] 1 0.08 0.06 0.04 -0.08 

B(K%7T°) 1.15 ±0.12 [32] 1 0.58 0.14 -0.95 

B(K^ri) 0.380 ± 0.060 [32] 1 0.10 -0.55 

B{K q s uj) 1.30 ± 0.30 [32] 1 -0.13 

B(K1tt°) 1.003 ± 0.083* [31] 1 



B(K~K + )/B(K-tt + ) 10.10 ±0.16 [28] 10.30 
bItt-tt + )/B(K-7t + ) 3.588 ± 0.057* [28] 1 



B. External Mixing Measurements 

1. Information on r 2 From Wrong-Sign D° — > K + tt~ and D° — > K^ + £~f£ 

Studies of wrong-sign D° — > K + n~~ from other experiments provide information on r 2 , y', 
and Rm- Without using D° decay times, these analyses are sensitive to the time- integrated 
wrong-sign rate ratio: 

Rws = r(/> _ (16) 

= r 2 + ry' + R u (17) 
= r 2 + v r 2 y cos 5 — v r 2 x sin 5 H — , (18) 

where Eq. (18) explicitly shows the dependence on the fit parameters r 2 , y, cos 5, x 2 , and 
x sin 5. By also measuring the decay times of the same D° candidates, these experiments can 
separate the three terms in Eq. (17) because each has a different proper time dependence: 



dN _ rf 
—— oc e 
dt 



r 2 + ry'Tt+ X ' 2+ A y, \ vt) 2 



(19) 



We use time-integrated i?ws measurements as a source of information on r 2 , which is 
otherwise poorly determined in our analysis, as shown in Table XV. Constraining r 2 leads 
directly to improved precision on cos S because it always appears as r cos 5 in Table I. Ex- 
tracting meaningful precision on r 2 from i?vvs requires information on x sin S and x 2 . When 
the external y and y' measurements described in Section IV B 2 are not included in the fit, 
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x sin 5 is poorly determined, so we take xsin<5 to be zero. For x 2 , which is also poorly de- 
termined in our analysis, we use Rm measurements from wrong-sign D° — > K^ + £^ui rates. 
Table XI shows the input i?ws and semileptonic Rm measurements, which are uncorrelated 
with all other measurements and are entered directly into the fit. 



TABLE XI: External mixing measurements used to constrain r 2 



Parameter Average (%) 
R ws 0.409 ± 0.022 [33-35] 

R M 0.0173 ± 0.0387 [36-39] 



2. Measurements of y and y' 



By combining external measurements of y and y' with the CLEO-c measurement of cos<5, 
the fitter can extract xsin<5 via xsinS = ycosS — y' . In doing so, it accounts for the fact 
that y' is determined simultaneously with r 2 and Rm, which are also functions of the fit 
parameters, and that cos 5 and its uncertainty are correlated with y (as discussed below in 
Section VI). 

External measurements of y come from two sources: comparison of D° decay times to 
CP eigenstates K + K~ and 7r + 7r~ with decay times to the flavor eigenstate K~ir + , and 
time-dependent Dalitz plot analysis of D — > KgTT + iv~, which is sensitive to both x and y. 
The latter x and y measurements are essentially uncorrelated, but x provides information 
on y when combined with measurements of Rm- 

For y', we use three sets of (CP- conserving) fit results from CLEO, Belle, and BABAR. 
Averages of these fits are shown in Table XII. The fit covariance matrices have been provided 
by the above collaborations and are also included in our analysis. 



TABLE XII: External measurements of y and y' with associated measurements of r 2 and x' 2 . 

Parameter Value (%) 

~y 0.662 ± 0.211 [6, 40-46] 

x 0.811 ± 0.334 [6, 46] 

r 2 0.339 ± 0.012 [4, 7, 47] 

y 1 0.34 ± 0.30 [4, 7, 47] 

x' 2 0.006 ± 0.018 [4, 7, 47] 



We note that, if we assume xsin<5 = 0, the y and y' measurements in Table XII by 
themselves provide an independent measurement of cos 5 = (y' + x sin S)/y = 0.52 ± 0.49, 
which is consistent with but less precise than the value found independently with CLEO-c 
data (see Section VI). 
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V. SYSTEMATIC UNCERTAINTIES 



All uncertainties except statistical uncertainties on the measured yields are considered 
to be systematic in origin. Systematic uncertainties on the fit inputs are included directly 
in the covariance matrix given to the fitter, which propagates them to the fit parameters. 
Uncorrelated yield uncertainties are discussed above in Section III. Below, we summarize 
the sources of correlated uncertainty and the measurements to which they apply. 

Final-state-dependent correlated systematic uncertainties, such as for tracking, particle 
identification (PID), and n° and K s finding efficiencies are the dominant uncertainties for 
the branching fractions, but, as discussed in Ref. [19], when no external measurements are 
included, they cancel in the DCS and mixing parameters. However, external measurements 
are, in general, not subject to the same correlated effects as our measured yields. Thus, 
while these external measurements reduce the statistical uncertainties on y and cos 5, they 
also introduce sensitivity to correlated uncertainties, and they increase the systematic un- 
certainties somewhat. 

In Table XIII, we list the correlated systematic uncertainties on reconstruction and PID 
efficiencies for final state particles, determined as Refs. [27, 52] describe. We also apply the 
efficiency corrections given in Refs. [27, 52] with the following modifications. For 7r° finding 
efficiency, we adjust the relative correction to —3.3% and inflate the systematic uncertainty 
to 4.0% to account for our inclusion of signal modes with more energetic 7r°s than in Ref. [27]. 
For f] finding efficiency, the relative correction appropriate for our selection criteria is —6.5% 
with an uncertainty of 4.0%. For K s finding efficiency, we include additional contributions 
of 0.3% and 0.4% for the flight significance and invariant mass requirements. We also include 
uncertainties on peaking background branching fractions, shown in Table IX. The signal 
and background efficiencies listed in this paper include the corrections discussed above. 

Table XIV shows mode-dependent uncertainties. Most are assessed by relaxing selection 
requirements and noting the resultant change in efficiency-corrected yield. We compute each 
of these uncertainties as the quadrature sum of the yield shift and its statistical uncertainty, 
where the latter quantity accounts for the large correlation between the two yields. 

The effect of modeling of the K^^it resonant substructure on the efficiency is determined 
by reweighting the simulated events according to the substructure observed in data and 
recomputing the efficiency. We assign an uncertainty for the modeling of multiple electron 
candidates by taking the difference between multiple candidate rates in simulated events 
and data. The efficiency of the extra track, n°, and r\ vetos in K^n modes is studied 
by comparing the frequency of these extra particles between simulated events and data. 
For semileptonic decays, we extrapolate the observed electron momentum spectrum below 
200 MeV/c based on the simulation. By sampling different form factor models used in the 
simulation, we estimate a relative uncertainty of 25% on this extrapolation. 

In our simulation, final state radiation (FSR) is generated by PHOTOS [53]. We measure 
efficiencies with and without FSR generation and assign 30% of the difference as a systematic 
uncertainty [27]. We also include a contribution corresponding to the change in efficiency 
when PHOTOS does and does not model interference among the radiated photons. 

Unlike FSR, initial state radiation (ISR) coherently shifts both M values upward in a 
signal DT. So, mismodeling of ISR would affect ST and DT efficiencies by the same fraction. 
Based on a study of excluding events in the ISR region (M > 1.87 GeV/c 2 ), we assign a 
conservative 0.5% uncertainty to all yields. 
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TABLE XIII: Correlated, fractional efficiency systematic uncertainties and the schemes for applying 
them in the data branching fraction fit. 



Source 


Uncertainty (%) 


Scheme 


Track finding 


0.3 


per track 


if ± hadronic interactions 


0.6 


per K^ 


Kg finding 


1.9 


per K s 


7T° finding 


4.0 


per 7T° 


r\ finding 


4.0 


per r\ 


dE/dx and RICH 


0.3 


per 7r ± PID cut 


dE/dx and RICH 


0.3 


per K± PID cut 


EID 


1.0 


per 



TABLE XIV: Correlated, mode-dependent fractional systematic uncertainties in percent for STs. 
An asterisk (*) marks those uncertainties that are correlated among modes. The schemes by which 
these uncertainties are applied to DTs are also given, along with the fractional DT uncertainty, 
Adt, on mode {^4, B} for ST uncertainties of a on mode A and (3 on mode B. 





AE 


ISR* 


FSR* Lepton Veto* 


Other 






0.5 


0.5 


1.2 


0.5 






K+K~ 


0.9 


0.5 


0.8 


0.4 


0.5 


K^ cos 9 cut 




1.9 


0.5 


1.7 


3.2 






Kg7T°TT° 


2.6 


0.5 






1.5 
0.7 


Kg daughter PID 
resonant substructure 


Kgir 


0.9 


0.5 










K s n 


5.5 


0.5 






0.3 
0.7 


rj mass cut 
B ( V - 77 ) [28] 


K s u 


1.2 


0.5 


0.8 




1.4 
0.8 


oj mass cut 

B(u -► tt+tt-tt ) [28] 


Xev e 




0.5 


0.3 




2.0 
0.7 


spectrum extrapolation 
multiple candidates 






0.5 






0.7 
0.3 
1.4 
1.6 
0.5 


background subtraction 
extra track veto 
signal shape 
extra 7r° veto 
r] veto 


Scheme 


per D 


per yield 


per D 


per ST 


per D 




Adt 


yV + (3 2 {a + P)/2 a + (3 





V« 2 + P 2 





VI. FIT RESULTS 

We tested the analysis technique with the full simulated sample of uncorrelated D°D° 
decays that was filtered to leave a 40% subset of events that mimic the effect of quantum 
correlations. The effective luminosity of this subset is 15 times that of the data sample. 
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The fit results showed satisfactory agreement with the input values, taking into account the 
statistical correlation between the signal efficiencies (measured in the full simulated sample) 
and the signal yields (measured in the 40% subset). 

In Table XV, we first show the results of a data fit that includes external branching 
fraction measurements from Table X, but not the external -Rws an d Ryi measurements in 
Table XI. The corresponding correlation matrix for the fit parameters is given in Table XVI. 
Because r 2 has a large uncertainty and a negative central value, we cannot extract cos 5; 
instead, we quote rcos5. Also, we fit for R M instead of x 2 . 

TABLE XV: Results from the fit with external inputs from Table X, but not Tables XI or XII. 
Uncertainties are statistical and systematic, respectively. Charge-averaged D° branching fractions 
are denoted by final state. 



Parameter 


Fitted Value 


N (io 6 ) 


1.062 ±0.024 ±0.011 


y (10- 3 ) 


-52 ±60 ±17 


r 2 (IO" 3 ) 


-24 ±16 ±12 


rcos<5 


0.089 ± 0.036 ± 0.009 


RM (IO' 3 ) 


2.0 ±1.2 ±1.2 


rxsind (10 -3 ) 


(fixed) 


B(K-tt+) (%) 


3.81 ±0.05 ±0.06 


B{K-K+) (IO- 3 ) 


3.86 ±0.06 ±0.06 


£(«+) (IO" 3 ) 


1.37 ±0.02 ±0.03 


B(K° S TT°TT ) (IO" 3 ) 


8.27 ±0.45 ±0.41 


B(K° s tt°) (%) 


1.13 ±0.03 ±0.03 


B(K° s r)) (IO" 3 ) 


4.36 ±0.15 ±0.27 


B(K° s u) (%) 


1.10 ±0.04 ±0.05 


B(X~e+u e ) (%) 


6.35 ±0.19 ±0.18 


B{Kl^) (%) 


1.02 ±0.03 ±0.02 


Xfi t / ndof 


26.1/44 



In order to control the uncertainty on cos 5, we include in our standard fit both external 
branching fractions from Table X as well as external -Rws an d Rm measurements from 
Table XI. In this fit, shown in Table XVII, we obtain a first measurement of cos 5, consistent 
with being at the boundary of the physical region and with a precision that is dominated 
by the CLEO-c {Kir, S±} yield measurements. Our value of y is consistent with the world 
average of 0.00662 ± 0.00211 (see Table XII). For this standard fit, we assume xsin<5 = 0, 
and the associated systematic uncertainty is ±0.03 for cos 5 and negligible for all other 
parameters. The correlation matrix for this standard fit is shown in Table XVIII. Our 
branching fraction results do not supersede other CLEO-c measurements. 

The likelihood curve for cos 5 is obtained by repeating this fit at fixed values of cos 5 and 
recording the change in minimum x 2 . We then compute C = e~^ x _x mm)/ 2 ; which is shown 
in Fig. 6a. It is highly non-Gaussian, so we assign asymmetric uncertainties by finding the 
values of cos<5 where A% 2 = 1 to obtain cos<5 = 1.03^^7 ± 0.06. For values of |cos5| < 1, 
we also compute £ as a function of \S\, and we integrate these curves within the physical 
region to obtain 95% confidence level limits of cos 5 > 0.07 and \5\ < 75°. 
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TABLE XVI: Correlation coefficients (%) for the fit in Table XV using both statistical and sys- 
tematic uncertainties. 
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TABLE XVII: Results from the standard fit (with inputs from Tables X and XI) and the extended fit 
(with inputs from Tables X, XI, and XII). Uncertainties are statistical and systematic, respectively. 
Charge-averaged D° branching fractions are denoted by final state. 



Parameter 


Standard Fit 


Extended Fit 


N (10 6 ) 


1.042 ±0.021 ± 


0.010 


1.042 ± 0.021 


i 


0.010 


y (io- 3 ) 


-45 ± 59 ± 


15 


6.5 ±0.2 


± 


2.1 


r 2 (10^ 3 ) 


8.0 ± 6.8 ± 


1.9 


3.44 ±0.01 


± 


0.09 


cos 5 


1.03 ± 0.19 ± 


0.06 


1.10 ±0.35 


± 


0.07 


X 2 (IO" 3 ) 


-1.5 ± 3.6 ± 


4.2 


0.06 ± 0.01 


± 


0.05 


xsin<5 (10 -3 ) 


(fixed) 




4.4 ±2.4 


± 


2.9 


K-TT+ (%) 


3.78 ± 0.05 ± 


0.05 


3.78 ± 0.05 


± 


0.05 


K'K+ (IO" 3 ) 


3.87 ± 0.06 ± 


0.06 


3.88 ± 0.06 


± 


0.06 


7T-7T+ (IO" 3 ) 


1.36 ± 0.02 ± 


0.03 


1.36 ±0.02 


± 


0.03 


K° s TT°ir (IO" 3 ) 


8.34 ± 0.45 ± 


0.42 


8.35 ± 0.44 


± 


0.42 


K S 7T° (%) 


1.14 ± 0.03 ± 


0.03 


1.14 ±0.03 


± 


0.03 


K° s r] (IO" 3 ) 


4.42 ± 0.15 ± 


0.28 


4.42 ±0.15 


± 


0.28 


K° s u (%) 


1.12 ± 0.04 ± 


0.05 


1.12 ±0.04 


± 


0.05 


X-e+u e (%) 


6.54 ± 0.17 ± 


0.17 


6.59 ±0.16 


± 


0.16 


Kl^ (%) 


1.01 ±0.03 ± 


0.02 


1.01 ±0.03 


± 


0.02 


Xfit/ndof 


30.1/46 




55.3/57 





The asymmetric uncertainties on cos 5 quoted above still do not fully capture the non- 
linearity of the likelihood. This non-linearity stems from the use of r 2 to convert r cos S into 
cos S, which causes the uncertainty on cos 5 to scale roughly like 1/r. Because r 2 is obtained 
from i? ws via Eq. (18), an upward shift in y lowers the derived value of r 2 (for positive 
rcos5). As a result, the uncertainty on cos 5 increases with more positive y, as shown in 
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TABLE XVIII: Correlation coefficients (%) for the standard fit in Table XVII using both statistical 



and systematic uncertainties. 
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TABLE XIX: Correlation coefficients (%) for the extended fit in Table XVII using both statistical 
and systematic uncertainties. 
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Fig. 6b. A second, smaller effect is that the sign of the correlation between r 2 and rcos5 is 
given by the sign of y. When y is positive, r 2 and rcos5 are anti-correlated, which tends to 
inflate the uncertainty on cos 5. 

We also perform an extended fit that determines xsin5 by including previous measure- 
ments of y and y' . Table XVII shows the results of this fit, which incorporates all external 
measurements, from Tables X, XI, and XII. The correlation matrix for the extended fit is 
given in Table XIX. The resultant value of y includes the CLEO-c measurement from the 
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3950108-008 




-0.5 0.5 1.5 0.0 0.5 1.0 1.5 2.0 

cos8 

FIG. 6: Standard fit likelihood including both statistical and systematic uncertainties for cos<5 
(a), and simultaneous likelihood for cos<5 and y (b) shown as contours in increments of la, where 
a = y Ax 2 . The hatched region contains 95% of the area in the physical region. 

standard fit, but the precision is dominated by the external y measurements. The overall 
uncertainty on cos 5 increases to ±0.36 because of the i?ws measurement, as discussed above. 
However, unlike the standard fit, the likelihood for cos 5 shown in Fig. 7a is nearly Gaussian. 
For x sin 5, we assign asymmetric uncertainties resulting in x sin 5 = (4.4i 2 ;£±2.9)xlCr 3 . By 
repeating the fit at various simultaneously fixed values of cos 5 and sin 5, we also determine 
5 = (221 

i2-ii)°- The corresponding 95% confidence level intervals within the physical region 
are cos<5 > 0.39, xsind G [0.002,0.014], and 5 G [-7°, +61°]. Performing this extended fit 
with y, x 2 , and x sin 5 fixed to zero results in a change in \ 2 of 25.1, or a significance of 5.0<r. 

In the standard fit, the large correlation between y and r 2 originates in the use of the 
i?ws measurement in Table XI to provide information on r 2 . Using Eq. (18), one sees that 
the precision on r 2 is driven by the precision on y. Hence, for positive rcos<5, r 2 and y are 
negatively correlated. Consequently, the cos 5 value obtained by dividing r cos 5 by Vr 2 is 
positively correlated with y. In the extended fit, these correlations are broken by the precise 
external measurements of y in Table XII, and the weak negative correlation between cos S 
and y results from Eqs. (10) and (11). 

We assess the fit's sensitivity to the assumed quantum correlation parameters in the 
peaking background subtractions by varying these parameters over their full allowed ranges. 
We find excursions of less than 1 x 10~ 4 in y and 0.003 in cos<5 for both the standard and 
extended fits. We assign a systematic uncertainty, included in Table XVII given by the size 
of these shifts. 
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40 

8 (degrees) 



FIG. 7: Extended fit likelihood including both statistical and systematic uncertainties for cos 5 
(a), :rsin<5 (b), 5 (c), and simultaneous likelihood for cos<5 and xs'm5 (d) shown as contours in 
increments of la, where a = \/ A% 2 . The hatched regions contain 95% of the area in the physical 
regions. For 5, the fit fails to converge beyond the limits of the plot. 



A. Analysis of Input Information 

In this Section, we probe the power of the individual fit inputs or groups of fit inputs 
by removing them from the standard fit one by one. We define the weight (or information) 
that each input contributes to fit parameter A to be the fractional change in \jo\ when it 
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is removed from the fit: 

j = V^(new) - Xjaj (standard) 
f ja\ (standard) 

Because of the aforementioned non-linearities, removing an input from the fit can result in 
a <j\ that is smaller than the standard value, leading to negative I\. 

Table XX shows I\ for A = y and cos<5. When I\ > 0, we also calculate the significance 
of the shift to be S\ = A A/ a/ cr^ (standard) — er^(new), where the minus sign accounts for 
the correlation between the two values of A. 

As expected from Eqs. (8), (9), and (12), information about y comes primarily from the 
{S±, e} DT yields, which supply 91% of the statistical power. In addition, the {Kit, e} yields 
provide the normalizing semileptonic branching fraction, and S± ST yields provide MBs ± - 
Because y is determined from ratios of these quantities, a meaningful measurement depends 
on the simultaneous presence of various combinations of inputs. For instance, removing the 
S± ST yields would also reduce the power of the {S±,e} DT yields because, then, Eq. (8) 
cannot be used. Table XX does not account for such double-counting of information; it 
simply shows I\ as each individual input is removed from the fit. For this reason, the sum 
of all I y exceeds 100%. 

This effect is even more pronounced with cos 5. As demonstrated in Eqs. (10) and (13), 
almost all the information about rcosS is contained in the {Kir, S±} DT yields, but only 
when they are combined with the £>s ± estimates provided by the S± ST yields. Furthermore, 
in order to obtain cos S from rcosS, we rely on the inputs from Table XI, which therefore 
have /cos s = 1. If either {Kit, S + } or {Kir, S-} yields are removed, then the asymmetries in 
Eqs. (10) and (13) cannot be formed. However, rcos5 can still be obtained with knowledge 
of MBk-tt+Bs ±) which comes from other inputs in the fit. 

A striking example of non-linearities in the fit occurs when {S'_,e ± } double tags are 
removed. Doing so causes y to fluctuate downward, which, in turn, lowers the uncertainty 
on cos 5, as discussed above. As a result, {S_, e^} double tags appear to have large negative 
information content. 

In performing this exercise, we also test for anomalous inputs that bias the fit results 
with undue weight. As Table XX shows, we find no evidence of instability from our choice 
of fit inputs. 



B. Purity of C-Odd Initial State 

The purity of the C-odd D°D° initial state may be diluted by a radiated photon, which 
would reverse the C eigenvalue of the D°D° system and thus bias the fit results. ISR, 
FSR, and bremsstrahlung photon emission are benign because they do not alter the relative 
angular momentum between the D° and D°, occuring either before the D°D° are formed 
or after they decay. One problematic process would be photon radiation from the ^(3770), 
resulting in a virtual + state that then decays to D°D°. This channel is suppressed, as 
there are no nearby + states available. Another possible contribution would be ^(3770) 
decay to a virtual D*°D°, which subsequently decays to _D°D°7. Theoretical estimates of 
these amplitudes indicate a ^(3770) — > D°/) 7 branching fraction of less than 10~ 8 [54]. 

We verify the absence of this effect in data by searching for same-CP DT signals, which 
would be maximally enhanced for the C-even configuration. As shown in Fig. 8 and Ta- 
ble XXI, all such modes have yields consistent with zero. To determine the C-even fraction 
of the D°D° sample, Mc+I-N'c-, we perform a variant of the standard and extended fits that 
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TABLE XX: Effect of removing inputs from the standard fit. For each input and fit parameter 
combination, we report the significance of the shift in central value S\ and the fractional information 

content I\. In cases where I\ is negative, we do not report the shift significance. 

Removed Input S y I y S cos s I C osS 
Single Tag Yields 
K-n+ +0.1(7 0% -O.Ocr 0% 
K+tt- +0.2(7 0% -O.Ocr 1% 

K + K~ 6% +0.4(7 41% 

it+tt- +0.6(7 1% +0.2(7 11% 

K°7r°7r° +0.5 2% +0.4(7 23% 

K° s 7r° +1.1(7 14% +0.5(7 72% 

K%r] -2.8(7 0% -31% 

ifg^ -2.5(7 2% — -46% 



Double Tag Yields 





+ 0.1(7 


1% 


+ 0.1(7 


8% 


K-TT + ,K + TT- 




-0% 




-2% 


KlT,S+ 




-0% 


-0.0(7 


48% 


Kir,S- 


-0.1(7 


0% 


+ 0.0(7 


38% 


S+, S- 


+0.9(7 


0% 


+0.2cr 


11% 


Kn,e ± 


-0.0(7 


11% 


-0.1(7 


1% 


S + ,e± 


-0.3(7 


37% 




-19% 


S-,e ± 


-0.9(7 54% 




-135% 


K°ir°,Kir 




-3% 


-1.3(7 


14% 






-0% +0.2a 


5% 


External or 


Other CLEO- 


c (*) Inputs 


Rws, Rm 


-0.6(7 


3% 




100% 


B(K-ir+) 


-1.2(7 


1% 




-12% 


B(K+K~) 


+ 1.4(7 


3% 


+0.2cr 


29% 


B{7t + 7T-) 




-0% 




-1% 


B(K° s ir ) 


+0.4(7 


1% 


-0.2cr 


3% 


B(K° s rj) 


+ 1.4(7 


0% 


+ 0.0(7 


5% 


B(K° s co) 
B(K° L ir )* 


+0.1(7 


0% 


-0.3cr 


0% 




-5% 


-1.2(7 


6% 



include the 15 same-CP DT yields and efficiencies shown in Table XXI. We then express each 
yield as a sum of C-odd and C-even contributions: JVjj = Afc-AjA^jRc- + Afc+AfA^Rc+, 
where Rc- and Rc+ are the functions of y, r, cos 5, and x sin S given in Table XXII. The 
results of these fits are shown in Table XXIII. The C-even fraction is found to be consistent 
with zero, with an uncertainty of 2.4%, and shifts of the fit parameters from the nominal fit 
results are negligible compared to the systematic uncertainties already assigned. 
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FIG. 8: Combined M distribution for fully reconstructed same-CP and {SLjS 1 -} DT 

modes and combined M^ iss distribution for {S' + , K®7r°} modes. Data are shown as points with 
error bars. For {<!?+, and {S^S-}, the shaded histogram shows events outside the signal 
region. For {S' + , K^ir }, the shaded histogram represents simulations of the peaking backgrounds 
D° -» ir°ir°, K%tt°, 7]ir°, and K*°ir°. 



VII. SUMMARY 



We employ a double tagging technique with quantum-correlated D°D° decays at the 
^(3770) to perform a first measurement of cos 5 = 1.031qi7 ^ 0.06, where the uncertainties 
are statistical and systematic, respectively. Within the physical region, we find \S\ < 75° at 
the 95% confidence level. By including external inputs on y and y' in the fit, we find an 
alternative value of cos 5 = 1.10 ± 0.35 ± 0.07, as well as x sin 
6 = (22^)°. 
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TABLE XXI: Same-CP double tag yields and efficiencies included in the C-even-allowed fit. 



Mode Yield Efficiency (%) 



K + K~,K+K- 


-2.: 


L± 1.5 ±3.1 


32.4 ± 1.7 


K + K-,ir + 7r- 


0.9 


± 


1.6 


± 1.1 


38.9 ± 2.2 


K + K~ ,K° s it -k 


0.0 


± 


1.0 


±0.0 


6.6 ±0.5 


7T + 7r~, 7T + 7T~ 


1.2 


± 


2.8 


± 4.8 


44.9 ±5.1 


7T+7T ,^7T U 7T U 


1.0 


± 


1.0 


±0.0 


9.2 ± 1.0 


Kg7r°7r ,Kg7r 7T 


0.0 


± 


1.0 


± 0.0 


1.5 ±0.2 




0.0 




1 


+ 00 


8.9 ±0.5 


t^O t^CU 


n n 


± 


1.0 


±0.0 




K° s Tr ,K° s LU 


1.0 


± 


1.0 


±0.0 


3.3 ±0.2 


K° sV ,K° sV 


0.0 


± 


1.0 


±0.0 


0.7 ±0.4 


K° sV ,K s u 


0.0 


± 


1.0 


±0.0 


1.0 ±0.2 


K®lu, K° s u 


1.0 


± 


1.0 


±0.0 


1.5 ±0.2 




1.7 


± 


2.1 


±0.0 


28.2 ±0.1 


K° L -K°,TT + -K- 


4.3 


± 


2.4 


±0.0 


38.5 ±0.1 


K° L 7T ,K s 7T 7r° 


2.7 


± 


2.3 


±0.0 


6.4 ±0.1 



TABLE XXII: Comparison of C-odd and C-even yields, normalized by N 'Af for ST modes i and 
by AfAjAj for DT modes to leading order in x, y, and r 2 , with y' = ycos5 — x sin 5 and 

y = ycosS + xsin<5. Charge conjugate modes are implied. 



Mode 


Rc~ 




Rc+ 


K-7T+ 


1 + 2ry cos 5 + r 2 




1 + 2ry cos 5 + r 2 


s+ 


2(1 -y) 




2(1 -y) 


S- 


2(1 + y) 




2(1 ±y) 


K-TT + ,K-7T + 


Rm 




2r(r + y') 
■ 2r 2 (l - 2 cos 2 5) + Ary 


K-TT + ,K + 7T- 


l + 2r 2 (l-2cos 2 5) 1- 


K-7T + ,S+ 


1 + 2r cos 5 + r 2 


[1- 


-2rcos5 + r 2 ] (1 - 2y) 


K~Tr + ,S- 


1 — 2r cos 5 + r 2 


[l + 2rcos5 + r 2 ] (1 + 2y) 


K~tt + , e~ 


1 




1 + 2ry 


S+, S+ 







2(1 - 2y) 


S-,S- 







2(1 + 2y) 


s + ,s^ 


4 







S+,e~ 


1 




l-2y 


S-, e~ 


1 




l±2y 



APPENDIX A: CORRECTED YIELDS 

For use by future experiments, we provide our efficiency-corrected, background-subtracted 
yields in Tables XXIV and XXV. The correlation matrix, including statistical and systematic 
uncertainties, for these yields appears in Tables XXVI-XXXI. 
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TABLE XXIII: Results from the fits allowing for a C-even component in the initial state, with 
inputs from Table XXI. Uncertainties are statistical and systematic, respectively. We also give the 



s in the fit results 


with respect to the nomina 


fits. 






Parameter 


Standard, C-even Allowed 


QVi; fl- 
oruit 


Extended, C—even Allowed 


oniit 


Mc- (10°) 


1 n A A i r\ non i n m T 

1.U44 ± 0.029 ± 0.017 


n nno 
±0.002 


1.1)44 ± (J.IMU ± U.U17 


n nno 
±0.002 


■Nc+I-Nc- l-W J 


— 1.0 ± ZO.U ± z.z 




1 A J- OQ I 1 O 

— 1.4 zt Zo.D ± i.y 




y (io ) 


—46 ± 59 ± 16 


n o 

—0.2 


o.o =b (J. 2 ± 2.1 


n n 

— u.u 


_2 /in— 3^ 
r (10 J 


8.0 ± o.o ± i.y 


n nno 
±0.008 


o.44 ± u.ui ± u.uy 


n nnn 
— U.UUU 


r 

COS 


1 no _i_ n 1 n _i_ n n£ 
1.02 ± u.iy ± 0.06 


n nno 
— 0.003 


1 in i n otr i n 

1.1U ± U.oo ± U.U/ 


— U.UUo 


™2 /in- 3\ 
X (1U J 


t C i A 1 1 O O 

— 1.5 ± 4.1 ± 3.8 


n o 

—0.2 


n n£? i n ni i n n c 

0.06 ± 0.01 ± 0.05 


i n n i 

±0.01 


xsmo (ID ) 


U (fixed) 




A A 1 O /I i o r\ 

4.4 ± 2.4 ± 2.9 


n n a 

— 0.04 


B(K 7T^J (7oJ 


o ^70 i r\ r\c i n nc 

3.78 ± 0.05 ± U.U5 


— 0.001 
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TABLE XXIV: Efficiency-corrected, background-subtracted single tag and double tag yields input 
to the fits. Uncertainties are statistical and systematic combined. 



Mode 




Yield 
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TABLE XXV: Efficiency-corrected, background-subtracted single tag and double tag yields input 
to the fits. Uncertainties are statistical and systematic combined. 
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Yield 



TT 1 TT 
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K^TT ,^ 
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K° S 7], K° s oj 
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X+e~u e 
X + e~u e 
X + e~v e 
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X + e~v e 
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K+ 



K 
K 

7T~7T + 

K° s tt°tt 
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Ve 
Ve 
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508 ± 105 



[17] Z. Z. Xing, Phys. Rev. D 55, 196 (1997) [arXiv:hep-ph/9606422]. 

[18] D. Atwood and A. A. Petrov, Phys. Rev. D 71, 054032 (2005) [arXiv:hep-ph/0207165]. 

[19] D. M. Asner and W. M. Sun, Phys. Rev. D 73, 034024 (2006) [Erratum-ibid. 77, 019901(E) 

(2008)] [arXiv:hep-ph/0507238]. 
[20] Q. He et al. [CLEO Collaboration], Phys. Rev. Lett. 95, 121801 (2005) [Erratum-ibid. 96, 

199903 (2006)] [arXiv:hep-ex/0504003]. 



32 



TABLE XXVI: Correlation coefficients (%), with statistical and systematic uncertainties, for the 
efficiency-corrected, background-subtracted yields in Tables XXIV and XXV. Coefficients of 100% 
are represented by a dash ( — ). 
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TABLE XXVII: Correlation coefficients (%), with statistical and systematic uncertainties, for the 
efficiency-corrected, background-subtracted yields in Tables XXIV and XXV. 
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TABLE XXVIII: Correlation coefficients (%), with statistical and systematic uncertainties, for the 
efficiency-corrected, background-subtracted yields in Tables XXIV and XXV. 
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TABLE XXIX: Correlation coefficients (%), with statistical and systematic uncertainties, for the 
efficiency-corrected, background-subtracted yields in Tables XXIV and XXV. Coefficients of 100% 
are represented by a dash ( — ). 
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TABLE XXX: Correlation coefficients (%), with statistical and systematic uncertainties, for the 
efficiency-corrected, background-subtracted yields in Tables XXIV and XXV. 
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TABLE XXXI: Correlation coefficients (%), with statistical and systematic uncertainties, for the 
efficiency-corrected, background-subtracted yields in Tables XXIV and XXV. Coefficients of 100% 
are represented by a dash ( — ). 











o 














o 


























+ 


+ 














+ 




t= 






























1 








3 


i 


«! 


_i_ 


o 


o 




3 






i 




% 










1 




1 


OCT) 


oco 


OCn OCT) 


+ 


i 


1 


OtO OCn 


ceo 




+ 


i 




1 


o 


o 








S< 


^ 












s< 


s< 




s< 


s< 


s< 


S< 


s> 

1 


+ 


+ 




t= 

O<0 


OCT) 


c- 

oco 


3 
oco 




1^ 




1^ 






1^ 


la" 
















if 






+ 


it 1 




if* 
~( 


it* 




1 

CD 


1 

CD 


1 


1 

CD 


1 

CI) 


1 

CD 


CD 


+ 

CD 


+ 

CD 


+ 

CD 


+ 

CD 


+ 

CD 


+ 

CD 


+ 

CD 


o 


o 

t= 


O 


o 

t= 


o 

t= 


o 

Is 


o 

t= 


o 




+ 


+ 


+ 


+ 


+ 


+ 


+ 


1 


1 


1 


1 




1 


1 




























i*S 


*S 




















sc; 


sc. 


sc. 


sc; 


sc; 






X + e v e ,K 7T+ 




8 


6 


3 


9 


4 


7 


36 


8 


8 


4 


10 


4 


7 


7 


7 


-0 


-0 


-0 


1 


1 


3 


X + e~9 e , K~K + 






1 


1 


2 


1 


1 


8 


2 


2 


1 


2 


1 


1 


1 


1 


-0 


-0 


-0 











Vi. — — — -L 











2 


1 


1 


6 


1 


2 


1 


2 


1 


1 


1 


1 


-0 


-0 


-0 








1 


X + e V e , K^TT U TT U 










6 


1 


4 


3 


1 


1 


5 


7 


1 


4 


7 


7 


-0 


-0 


-0 


10 


2 


7 


v+ r/0 

X^e v e ,Kgir 












2 


7 


9 


2 


2 


9 


14 


2 


7 


12 


12 


—0 


— 


— 


17 


5 


12 


T/-4- j-^n 

X + e v e ,Kgr] 














1 


4 


1 


1 


1 


2 


5 


1 








-0 


-0 


-0 


1 


6 


1 


















7 


1 


2 


6 


8 


1 


7 


8 


8 


-0 


-0 


-0 


10 


3 


9 


X e+v e ,K+TT 


















8 


8 


4 


10 


4 


7 


6 


6 


-0 


-0 


-0 


1 





2 


X^e+v e ,K^K+ 




















2 


1 


2 


1 


1 


1 


1 


-0 


-0 


-0 











~\r -4- -4- 
X e^U ei 7T 7T^ 






















1 


2 


1 


2 


1 


1 


-0 


-0 


-0 








1 


,r _i_ n n 

X e + v e ,K%n u Tr u 
























10 


1 


G 


11 


11 


-0 


-0 


-0 


15 


4 


11 


■\r — _|_ T-j-fi n 


























2 


9 


13 


14 


-0 


-0 


-0 


20 


5 


14 


X e^v e ,KgT] 




























1 








-0 


-0 


— 


1 


5 


1 


X~e + v„ Kiu) 






























8 


8 


-0 


-0 


— o 


10 


3 


10 


































25 


-0 


-0 





26 


7 


19 


Ar£7r°,A' + 7r- 


































-0 


-0 





26 


7 


20 












































-0 


-0 


-0 











































-0 


-0 


-0 










































-0 


-0 


-0 












































9 


24 














































7 

















































38 



